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I. INTRODUCTION 

The two-dimensional layer of graphene is a most inter- 
esting object for the theoretical and experimental study 
of quantum phenomena. It has quite high mobility tun- 
able by a gate voltage and a unique band structure simi- 
lar to the energy spectrum of massless Dirac fermions ^^i^ 
Such a quasi-relativistic nature of free carriers affects 
the transport phenomena including the weak localization. 
The weak localization (antilocalization) is caused by the 
constructive (destructive) interference of electron waves 
traveling in opposite directions.^ The interference is sup- 
pressed by an external magnetic field, which leads to a 
magnetoresistance in classically weak fields. The anoma- 
lous magnetoresistance in graphene is the subject of in- 
tensive experimental^ — and theoretical^— study dur- 
ing the last few years. It was experimentally shown that 
increasing the carrier density and decreasing the tem- 
perature leads to a transition from weak antilocalization 
to weak localization in graphene Such a behavior is 
attributed to the trigonal warping of electron spectra in 
valleys, which suppresses antilocalization, and intervalley 
scattering, which restores localization.^^ 

Depending on the ratio between the magnetic length 
Ib and the mean free path i one distinguishes between 
two regimes of weak localization. In very low magnetic 
fields {Ib ^ the main contribution to the magne- 
toresistance comes from large diffusion- like trajectories of 
electrons with the size L ^ £. This is the diffusion regime 
of weak localization. With the field increase {£b ^ 
the role of large trajectories is suppressed and the weak- 
localization correction to conductivity is determined by 
electron trajectories with few scatterers. Such a non- 
diffusion regime is realized in high-mobility structures at 
rather small magnetic fields — To correctly extract ki- 
netic parameters of carriers, such as the phase breaking 
time, from the magnetoresistance measurements one has 
to analyze experimental data in the whole range of classi- 
cally weak magnetic fields. Previous calculations of mag- 
netoresistance in graphene were carried out only for the 
diffusion regime though experimental data in Ref.— in- 
dicate that the magnetic length may reach the mean free 
path in the magnetic field as small as 30 mT. The goal 
of this work is to develop the theory of weak localiza- 
tion beyond the diffusion regime and derive the quantum 



correction to conductivity. 

The developed theory can be also applied to two- 
dimensional systems formed at surfaces of bulk topologi- 
cal insulators, such as Bi2Se3, Bi2Te3, and Bi2Te2Se (see 
Refs. p2ll23j and references therein), and in HgTe quan- 
tum wells of critical width. The excitations in these 
systems are described by the effective Hamiltonian simi- 
lar to one for the free carriers in graphene. We note that 
the Dirac cones in graphene are situated at the K points 
of the two-dimensional Brillouin zone and degenerate in 
spin. In contrast, the similar energy spectrum of carriers 
in topological insulators is formed by spin-orbit interac- 
tion, with the Dirac point being situated at the T point of 
the two-dimensional Brillouin zone. Experimental data 
indicate that the magnetoresistance in Bi2Te3 (Ref. [25l ) 
and Bi2Si3 (Ref. 26) as well as in HgTe quantum well^ is 
positive (antilocalization) as expected for massless Dirac 
fermions. 



II. ORIGIN OF WEAK ANTILOCALIZATION 
IN GRAPHENE. 

Technically, the calculation of quantum corrections to 
conductivity involves the computation of Cooperons^^ 
In the case of multi- valley systems with negligible spin- 
orbit coupling, such as graphene or silicon, the Coore- 
pons are derived from an integral equation with the 
kernel (]/cK/3(k, k')Ky(5(— k, — k')), where the indices a, 
/3, 7, and S enumerate valleys, k and k' are the elec- 
tron wave vectors measured from the valleys centers, 
and the angle brackets denote averaging over the posi- 
tions of scatterers. We note that in a single- valley sys- 
temdue to the time inversion symmetry, one may take 
the orbital wave functions of free electrons in the form 
V^k = V^*k, which yields V(-k,-kO = F(k^k). Then, 
the correlator (l^(k, k') V(— k, — k')) equals to the correla- 
tor ^|V(k, k')|^^ determine the single-particle transport. 

It leads to a diffusion pole in the Cooperon equation and, 
therefore, to an enhancement of backscattering (weak lo- 
calization). In contrast, Vaa{—^^—^^) 7^ ^^^(k^k) gen- 
erally in multi- valley systems, there are no diffusion poles 
in intra- valley Cooperons, and weak localization is absent 
in the approximation of independent valleys However, 
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the special form of scattering amplitude may result in a 
diffusion pole even in the one-valley model. 

To consider this phenomenon for graphene we neglect 
the small trigonal warping of electron spectrum and take 
the effective Hamiltonian describing electron states in 
each valley in the form 



H 



Px- iPy 

Px + iPy 



(1) 



Here, v is electron velocity and p = (px^Py) is the mo- 
mentum operator. The wave functions of conduction 
band may be chosen in the form 



*k(r) 



ikr 



(2) 



where (p is the polar angle of the in-plane wave vector 
k. We consider n-doped graphene and assume that the 
valence-band states do not contribute to low-temperature 
conductivity. 

The intra- valley electron scattering from a short-range 
impurity {Ai symmetry) is described by the Hamiltonian 



SH 



1 
1 



(3) 



where vj is the impurity position. It follows from Eqs. (j2j) 
and ([3]) that the matrix element of scattering has the form 



V(k^k) =a 



i(k-k')Tj 



(4) 



^^i(^-^')/2cos[((^-(p')/2] . 



One can see that the direct back scattering from an impu- 
rity is suppressed and, what is more important for quan- 
tum effects, the scattering introduces the phase {(p—(p')/2 
to the electron wave function. Therefore, an electron 
traveling clockwise along a closed path and finally scat- 
tered back gains the additional phase 7r/2 while the elec- 
tron traveling in the opposite direction gains the phase 
— 7r/2. The phase shift of tt between these two waves 
results in a destructive interference and, hence, in the 
antilocalization of carriers. 

We note that other forms of scattering amplitude, e.g., 
y(k',k) oc a + bex.-p[i{(p — cp')] with a ^ lead to the 
phase gain which depends on the particular trajectory 
even for closed paths. Averaging over the trajectories 
destroys the wave interference and results in no quantum 
corrections to conductivity (see also Refs. 0, [H, andlii). 
Similar arguments apply for a surfaces of bulk topological 
insulators and HgTe QWs with the Dirac-like spectrum 
of carriers. 



III. CALCULATION OF 
MAGNETOCONDUCTIVITY. 

To calculate the interference corrections to conductiv- 
ity in the perpendicular magnetic field we use the dia- 



gram technique. The retarded and advanced Green func- 
tions of conduction electrons in one valley in graphene 
have the form 



eF-SN± ih/{2T) ± ih/{2T^) ' 



(5) 



Here, Sf = vfikp is the Fermi energy, kp is the Fermi 
wave vector, Sn = huJcVN ^ ujc = \/2v/^b-, ^ is the quan- 
tum relaxation time, 1/r = niO?kF / {2h^v)^ rii is the im- 
purity density, is the phase breaking time, we assume 
that T(j) ^ r, ^ N.kyij^) are the two-component wave func- 
tions given in the Landau gauge by 



'4'N,ky{T^) are the standard functions of a two-dimensional 
charge particle in magnetic field, and N and ky are the 
quantum numbers. Assuming that the magnetic field is 
classically weak, i.e., ujcT ^ 1, we obtain 



G^'^(r,r') = exp 



.ix + x')(y-y') 



2l\ 



G^>^(r-r'), (6) 



where G^'^{r — r') are the Green functions at zero field. 



^R,A,„^ _ exp[-p/(2f)±i(fcF/3 + 7r/4)] „^ 



y^27rp/kF hv 



i' = £/ (1 -^t/T(I))^ £ = 'ur is the mean free path, cr is 
the vector of Pauli matrices, n = p/p is the unit vector 
pointing along p, and it is assumed that kpp ^ 1. We 
note that the Green functions are matrices 2x2 owing 
to the matrix form of the Hamiltonian ([1]). 

The Cooperon C(r, r') may be presented as a matrix 
4 X 4 in the basis of direct product of the states (1,0)^ 
and (0, 1)^ inside the valley. It is found from the integral 
matrix equation 



C(r,r') = «;P(r,r') + j P{y,y^)C{y^,y') dn 



(7) 



with the kernel 



P(r, r') = wG^(r, r') ® G^(r, r') = 



Po(r,r') 



/ 1 



iri- —iri- n_ 



\ 



m+ 1 —1 —in 
m+ —1 1 m_ 
\ m+ — m+ 1 / 



(8) 



where w = niO? and 

-\v-v'\/l' 



Po(r,rO 



27r^|r-r 



-exp 



.{x^x'){y-y') 
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The matrix elements Pj,^(r, r') in Eq. (j8j) are obtained and L}^_m Laguerre polynomials. Having de- 

from the components G^^(r, r')G:^^(r, r') where the in- 



dices u and /i enumerate the states (a, 7) and re- 
spectively. 

A standard method to solve Cooperon equations is to 
expand the kernel in series of the wave functions ^at/c^ (r) 
of a particle with the charge 2e in the magnetic field 
To solve matrix Eq. ([7]) we modify the approach and in- 
troduce the basis matrices 



^N>i,ky{r) = 
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(9b) 
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One can see that J <l>jy ^^(r)<l>Ar',/e^(r)<^r 
with I being the unit matrix 4x4. 

Direct calculation shows that the matrix kernel P(r, r') 
is expanded in series of <l>Ar,fe^(r) as follows 



where 
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are coefficients given by the integrals 
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[x')x^dx, (12) 



composed the kernel P(r, r'), one can readily find that 
the solution of Eq. ([7]) for the Cooperon has the form 

C(r,r')= ^*iV,fc„(r)(I-Piv)-'Piv<fe„(r'). (13) 

N,ky 

The weak localization correction to conductivity 
has two contributions corresponding to the standard 
diagrams 

(14) 



cr = cr, 

The term cFa is given by 

7T 



CTa 



Tr 



P(r,rOC(^)(r^r)' 



drdr^ 



(15) 



where P(r, r') = Ja^(r, r')(g) Ja^(r, r'), J(r, r') is the current 
vertex, Ci^)(^r\r) = J P(r,ri)C(ri,rOdri, and the valley 
and spin degeneracy is taken into account. In graphene, 
the vertex can be presented in the form 



J(r,rO = 2ev J G^(r, n) tr G^(ri, r^^iri , 



(16) 



where the factor 2 stems from the difference between the 
quantum and transport relaxation times. Straightfor- 
ward calculation shows that the conductivity correction 
aa assumes the form 



a 2 /92 ^ 



(17) 



P(r, r') = J2 ^N,k, (r)P^*^,fc„ (r') , (10) where 



n 
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Similar procedure can be carried out to derive the con- 
ductivity correction (75 corresponding to nonbackscatter- 
ing interference effects. The calculation is more compli- 
cated, because it requires the expansion of each of the 
current vertices J(r,r') in the series of ^N,ky^ and yields 

^6 = Tr 1 Q^UQoil - Po)-'Po 



N=0 



where Qn are the matrices (here we assume P^ 
for M > N) 



(M) 
N 




Qn>i 





i- 
i 





p(i) 

N-1 
d(2) 











:)(2) 
N 



N 



V2 V2 

p(l) 



p(l) 

2 

p(3) 



(19a) 



4 




10-4 10-3 10-2 10-1 10° 10^ 10^ 10^ 10^ 



FIG. 1. Weak localization correction to conductivity as a 
function of magnetic field. Curves 1, 2, 3, and 4 are plotted 
for t/t0 = 10-^, 10-^, 10-^ and 10~^ , respectively. Dashed 
lines show the weak-localization corrections to conductivity 
at zero magnetic field 
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FIG. 2. Magnetic field dependences of the contributions a a 
and at. Curves are plotted for the same ratios r/r^ as in 
Fig. 1. Inset illustrates the dependences (j{B) — cr(0) on mag- 
netic field calculated after accurate Eqs. (pT|) and ([T8|) (solid 
curves) and in the diffusion approximation (dashes curves) 
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Equations ([11]), ([T7|), and (fT8|) describe the quantum 
corrections to conductivity in the whole range of classi- 
cally weak magnetic fields. 



IV. RESULTS AND DISCUSSION 



Figure 1 presents the magnetic field dependence of the 
conductivity correction (j{B) calculated for different ra- 
tios between the relaxation time r and phase breaking 
time Tcf). One can see that the conductivity correction is 
positive in the whole range of magnetic fields (weak an- 
tilocalization) and monotonously decreases with B giving 
rise to negative magnetoconductivity. Such a behavior is 
in accordance with the qualitative analysis of intravalley 
interference effects presented above. The magnetic field 
dependences of the contributions cia and are plotted in 
Fig. 2 illustrating that the contributions are of opposite 
sign and comparable in magnitude. 

Equations (fT7|) and (fT8|) allow us to analyze the behav- 
ior of conductivity in low and high magnetic fields. In the 
low- field limit, the contributions to magnetoconductivity 



(20) 



where F2{x) = lnx+'0(l/2+l/x), i/j is the digamma func- 
tion. These dependences are in agreement with the result 
of Ref. 112| obtained in the diffusion approximation. In 
the inset of Fig. 2, we compare the magnetoconducitivity 
(j{B) — cr(0) calculated after accurate Eqs. (p!7|) and (p!8|) 
(solid curves) with that calculated in the diffusion ap- 
proximation Eqs. (|2Q]) (dashes curves). One can see that 
the diffusion approximation describes the magnetocon- 
ductivity in low fields, where > ^- In higher magnetic 
fields, the diffusion approximation is not valid and one 
has to use the microscopic theory developed in this pa- 
per to describe the magnetoconductivity. Particularly, 
the ratio 2:1 between \cf a{B) — cf a{S))\ and |<J5(5) — (75(0)! 
obtained in the diffusion approximation does not hold 
any more. In the high- field limit, i.e., Ib <^ ^, we may 
keep only first order in Ib/^ terms in Eq. ([12]), which 
gives cr oc 1 / y/B. 

The absolute value of the weak-localization correction 
to conductivity in zero magnetic field can be obtained 
from Eqs. (fT7j) and (fTSj) by taking a formal limit Ib ^ 
0. The calculated cr(0), 0-^(0), and 0-5(0) are shown in 
Figs. 1 and 2 by dashed hues. The correction (j(0) is 
determined by the phase breaking time, which can be 
used to experimentally study the dependence of upon 
temperature or carrier density. 

Finally, we note that the presented theory can be ap- 
plied to describe the weak localization of two-dimensional 
Dirac fermions in conducting surface of bulk topological 
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insulators. In such systems, there is one two-dimensional 
channel and the conductivity corrections are given by 
Eqs. (dll) and divided by 4. 

To summarize, we have developed the theory of weak 
localization for graphene beyond the diffusion regime and 
calculated the interference corrections to conductivity in 
the whole range of classically weak magnetic fields. The 
theory will allow one to better describe experimental data 
and more precisely determine the phase relaxation time 
of electrons in graphene. 
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